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ABSTRACT 

We report the discovery of two new halo velocity groups (Cancer groups A and 
B) traced by 8 distant RR Lyrae stars and observed by the Palomar Transient 
Factory (PTF) survey at RA ~ 129°, Dec ~ 20° (1 ~ 205°, b ~ 32°). Located 
at 92 kpc from the Galactic center (86 kpc from the Sun), these are some of 
the most distant substructures in the Galactic halo known to date. Follow-up 
spectroscopic observations with the Palomar Observatory 5.1-m Hale telescope 
and W. M. Keck Observatory 10-m Keck I telescope indicate that the two groups 
are moving away from the Galaxy at v^ sr = 78.0 ± 5.6 km s _1 (Cancer group 
A) and Vg Sr = 16.3 ±7.1 km s" 1 (Cancer group B). The groups have velocity 
dispersions of = 12.4 ± 5.0 km s _1 and = 14.9 ± 6.2 km s~ , and 

are spatially extended (about several kpc) making it very unlikely that they 
are bound systems, and are more likely to be debris of tidally disrupted dwarf 
galaxies or globular clusters. Both groups are metal-poor (median metallicities 
of [Fe/H] A = — 1.6 dex and [Fe/H] B = — 2.1 dex), and have a somewhat uncertain 
(due to small sample size) metallicity dispersion of ~ 0.4 dex, suggesting dwarf 
galaxies as progenitors. Two additional RR Lyrae stars with velocities consistent 
with those of the Cancer groups have been observed ~ 25° east, suggesting 
possible extension of the groups in that direction. 
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Subject headings: stars: variables: RR Lyrae — Galaxy: halo — Galaxy: kine- 
matics and dynamics — Galaxy: structure 



Introduction 



In the last de cade, deep and w ide-area astronomical surveys, such as the Sloan Digital 
Sky Survey fSDSS: lYork et alihoooh and Two-Micron All Sky Survey, (2MASS: lskrutskie et al 
20061 ) have enabled detection of substructures in the halo both as stellar overdensities in 
space and as mov ing groups. The tidal streams of the disrupting Sagittarius dwarf galaxy 
fllbata et al.lll994f) are the best exampl e of such substructures, with streams wrapping around 



density ("Clump J" in 



Sesar et al. 



most of the sky (IMaiewski et al.ll2003f). Other known substructures include the Pisces Over 



2007 ). and several other overdensities and streams^. 



2007: Watkins et al. 2009; Sesar et al. 2010a). the GD-1 



strea m ( Grillmair fc Dionatos 2006 ). the Orphan stream ( GrillmahT2006 : Belokurov et al. 



While many halo substructures have been discovered so far, the motivation for find- 
ing and characterizing more of them remains strong. For example, at larger galactocentric 
distances (-Rgc ^ 15 kpc), the Galactic gravitational potential is dominated by dark mat- 
ter and observations of more distant tidal streams can be used to constrain the shape 



orientation, and ma ss of the Milky Way's dark matter halo (e.g., lLaw fc Majewskil 12010 



Watkins et al.ll2010l ). Other motivations for completing the census of halo sub structures are 



to better constrain the severity of the so-called "missing satellite problem" (jKlypin et al. 



19991 : iMoore et al.lll999l ) and the contribution of the accretion of dwarf satellite galaxies to 
the formation of the Galactic halo. 

Finding distant halo substructures, however, is not an easy task. One approach to 
this problem is to look for spatial groups of (non-variable) sources that are consistent in 
color-magnitude space with, for example, an old population of stars at a fixed distance (e.g., 
Willman et al.ll2005t lGrillmairll2006t iBelokurov et al.l 120071 ). While this technique has been 
employed quite successfully in the past, its strongest limitation when looking for distant 
halo substructures using ground-based imaging data is the separation of stars from galaxies 
at faint magnitudes. At magnitudes fainter than r ~ 21 (corresponding to a heliocentric 
distance of ~ 30 kpc for metal-poor main-sequence turn-off stars), the number of field Milky 
Way stars per unit magnitude decreases with increasing magnitude (i.e., the number den- 
sity profile of Galactic halo stars steepens with distance from the Galactic center from a 



A fairly extensive list of halo streams and overdensities can be found at 



http : //homepages . rpi . edu/~newbeh/mwstructure/MilkyWaySpheroidSubstructure . html . 
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powe r-law with index n = —2.6 to that with n < —3.8; iSesar et al.l 120111 ; iDeason et al. 
201 ll ). while the number of galaxies per unit magnitude increases. In addition, morpho- 
logical separation o f stars and galaxies (e.g., SDSS resolved-unresolved source classification; 
Lupton et al.l 120021 ) becomes increasingly unreliable towards faint magnitudes as the signal- 
to-noise ratio decreases. The increasing contribution of galaxies in point-source catalogs at 
faint magnitudes causes a sharp increase in the noise from which the signal of a true low 
surface brightness stellar system must be detected. 

Instead of looking for clustering in samples of non- variable tracers, which may get con- 
taminated by galaxies at faint magnitudes, we should look for clustering in samples of variable 
tracers with distinctive light curves (e.g., RR Lyrae stars) that are difficult to confuse with 
galaxies. RR Lyrae stars have several advantages when used to map the Gal actic halo be - 
cause they i) are old, evolved stars and therefore trace old stellar populations ( jSmithlll995l ); 
ii) have distinct light curves which make them easy to identify given multi-epoch observa- 
tions (peak-to-peak amplitudes of V ~ 1 mag and periods of ~ 0.6 days); and iii) are bright, 
standard candles (My = 0.6 mag at [Fe/H] = —1.5 dex, with ~ 7% uncertainty in distance) 
that can be detected at large distances (5-120 kpc for 14 < V < 21). The steepening of the 
stellar density profile beyond 30 kpc is actually beneficial for searches that employ RR Lyrae 
stars because it reduces the pool of stars that can form false spatial groups, and therefore 
increases the contrast between stars associated with halo substructures and stars associated 



with the smooth halo (e.g., the Pisces Overdensity shown in Figure 11 of lSesar et al.ll2010al ). 
Due to the steepening of the stellar density profile beyond 30 kpc, groups of distant RR 
Lyrae stars are more likely to be real halo substructures than chance associations of stars. 
While the probability of chance association at large distances and small angular scales is 
very small, it is not zero and spectroscopic followup is still needed to test whether the stars 
in a spatial group also form a moving group. 

In this paper we use t he line of reasoning presented above, and successfully used by 
Kollmeier et al.l (120091 ) and ISesar et al.l (|2010b|), as a motivation to study a group of 8 RR 
Lyrae stars found during a preliminary search for halo substructures in regions of the sky 
observed by the Palomar Transient Factory survey^. The data set and the algorithm used 
to select RR Lyrae stars are described in Sections [2] and [3J The spectroscopic observations, 
their reduction, and measurement of line-of-sight velocities and metallicities are described 
in Section [5j The results are discussed in Section [6] and our conclusions are presented in 
Section [71 



http : //www . astro . caltech . edu/ptf 
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2. Overview of the Palomar Transient Factory Survey 



The Palomar Transient Factory (PTF; lLaw et al.ll2009l . iRau et al.ll2009l ) is a synoptic 
survey designed to explore the transient sky. The project utilizes the 48-inch Samuel Oschin 
Schmidt Telescope on Mount Palomar. The telescope has a digital camera eq uipped with 
12 CCDs (one of which is not active), each 2Kx4K pixels ( IRahmer et al.ll2Q08h . Each PTF 
image covers 7.26 deg 2 with a pixel scale of 1.01". By the end of 2011, the PTF observed 
~ 8000 deg 2 of sky at least 30 times in the Mould- R filtejfl (hereafter,the i?-band filter) and 
~ 1800 deg 2 in the SDSS g' filter. 

Some regions of the sky, known as the "high-cadence" regions in PTF, have been ob- 
served more than 100 times. One such region on which we focus in this paper is a 45 deg 2 
field roughly centered on the Praesepe3 cluster (08 h 40' 24", +19° 41'), which has bee n 



imag ed more than 200 times by PTF in the i?-band (lAgueros et al.l 1201 ll ; iPolishook et al. 



20121 ). Relative to other, less observed areas, the Praesepe region is interesting because, 
due to chance and a large number of observations, it has a non-negligible number of images 
that are deeper than the nominal 5a faint limit at R = 20.6 (i.e., images that were taken 
in better than average conditions). In turn, these deeper images increase the likelihood of 
finding distant RR Lyrae stars and halo substructures, if any are present. 

We use PTF data processed by the photometric pipeline, which is hosted by the Infrared 
Processing and Analysis Center (IPAC). This pipeline perf orms final image reduction, source 



extra ction, and photometric and astrometric calibration (jGrillmair et al.l I201CH : lOfek et al. 



2012c Laher, R. et al. 2012, in prep). The photometric uncertainty provided by this pipeline 
is smaller than ~ 0.01 mag for R < 16 sources and increases to 0.2 mag at R = 20. 6. 



he a lgorithm used for phot ometric cali b ration is b ased on the one prese nted in iHoneycutt 



Ofek et al. 


( 


2011 


) and 


Levitan et al. 


f 


2011) 



reference UCAC-3 astrometric catalog (IZacharias et al.ll2010l ). the astrometric precision of 
PTF coordinates is about 0.1" in right ascension and declination. To enable fast processing 
and manipulat ion, the da ta used in this paper were stored in the Large Survey Database 
format (LSDfl Ijuridhoilh . 



3 The PTF Mould- R filter is similar in shape to the SDSS r-band filter, but shifted 27 A redward. 
4 Also known as the Beehive Cluster and M44. 
fottp : //lsddb . org| 
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3. Selection of RR Lyrae Stars 

In Figured] we show the spatial distribution of 8 distant (galactocentric distance Rqc > 
84 kpc) a6-type RR Lyrae (RRab) stars observed by PTF in the Praesepe region. These 
stars were selected during a preliminary search for RR Lyrae stars in the av ailable PTF 
observa tions. The selection procedure we used is similar to the one outlined by ISesar et al. 
(]2010al ) and is briefly described below. The completeness of this sample is estimated in 
Section HI 

First, we searched for all PTF sources positionally associ ated within 1.5" of point-like 
sources in the SDSS Data Release 8 catalog (lAihara et al.ll201ll ) that ha ve single-epoch S PSS 
colors consistent with colors of RR Lyrae stars (Equations 6 to 9 from lSesar et al.ll2010al ): 



0.75 <u- g < 1.45 
-0.25 < g-r < 0.4 
-0.2 < r-i < 0.2 
-0.3 < % - z < 0.3. 



(1) 
(2) 
(3) 
(4) 



In the above eq uations, SPSS colors have been corrected for ISM extinction using the dust 
map created by ISchlegel et al.l (119981 ) . 



Variable sources were selected from this color-selected sample by requiring that PTF 
i?-band light curves have root-mean-square (rms) scatter greater than 0.1 mag and x 2 P er 
degree of freedom greater than 3 (calculated with respect to the mean magnitude associ- 
ated with a given source). The five most likely periods of variability were obtained for each 
variabl e PTF source by running an imp lementation of the Supersmoother period-finding algo- 
rithm (lFriedmanlll984l ; lReimannlll994j ) on the source's i?-band light curve. If Supersm oother 
retur ned one or more periods in the 0.2-0.9 day range (typical of RR Lyrae stars; ISmith 
19951 ). then the variable source's light curve was phased (period- folded) with each p eriod 
and SPSS r-band RR Lyrae light curve templates constructed by ISesar et al.l (j2010al ) were 
fitted to phased data. 

The light-curve fitting was performed to minimize the robust goodness-of-fit cost func- 
tion defined in Equation [5] in the least-square sense, with the epoch of the peak brightness 
0o, peak-to-peak amplitude A R , and peak brightness R as free parameters. The quality of 
a template fit was defined with a x 2 -like parameter 



£ — median ( | Tn observed ra te mpiate\/ ^observed) ■> 



(5) 



where m b serve d and e b se rved are the observed R band magnitude and its uncertainty, m te mpiate 
is the magnitude predicted by the template, and i = l,N obs , where N obs is the number of 
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observations. Here we used the median to minimize the bias in ( due to rare observations with 
unreliable (non-Gaussian) errors (e.g., due to image artifacts, cosmic rays). The template 
with the lowest ( value was selected as the best fit, and the best-fit parameters were stored. 
Finally, we visually inspected template-fitted light curves to tag sources as RR Lyrae stars. 

The positions and light curve parameters of RR Lyrae stars observed in this work 
are listed in Table [1] and examples of phased i?-band light curves are shown in Figure |2j 
The periods and epochs of maximum brightness are uncertain to within a few seconds and 
minutes, respectively, allowing the phase of pulsation to be determined with better than 1% 
uncertainty. An accurate estimate of phase is important when subtracting the velocity due 
to pulsations from the measured radial velocity, as we explain in Section 15.31 In addition 
to 8 RR Lyrae stars discussed in this work (RR1 to RR8), Table [T] contains two RR Lyrae 
stars that were spectroscopically followed-up after initial submission of this paper (RR9 and 
RR10). These two stars were not used in the analysis presented below. However, their 
position and velocities are mentioned in Section [7J as these stars may indicate an eastward 
extension of the Cancer groups. 

The heliocentric distances of RR Lyrae stars, d, are estimated from measurements of 
their flux-averaged and extinction-corrected i?-band magnitudes, (R) , and by assuming that 
Mr m My. The Mr My approximation is uncertain at ~ 0.04 mag level, as estimated 
from the r ms scatter of the diff erence of flux-averaged r and V band magnitudes of RR Lyrae 
stars from ISesar et al.l ( )2010al ) (see their Table 3) . The abs o lute m agnitudes of RRab stars 



in the Johnson l/-band were calculated using the IChaboyerl (119991 ) My — [Fe/H] relation 



M v = (0.23 ± 0M)[Fe/H] + (0.93 ± 0.12), (6) 

where [Fe/H] is the spectroscopic metallicity measured in Section [5~4l The fractional error 
in the heliocentric distance due to uncertainty in Equation El spectroscopic metallicities, and 
RR Lyrae evolution off the zero-age horizontal branch is estimated at 7% (~ 6 kpc error in 
distance at 86 kpc). The apparent (R) magnitudes of the observed stars range from 19.9 
mag to 20.3 mag, corresponding to heliocentric distances of 77 kpc and 96 kpc, respectively. 



4. Completeness of the RR Lyrae Sample 

The selection procedure described in the previous section yielded 8 RRab stars at dis- 
tances from 77 to 96 kpc from the Sun in the 45 deg 2 region of the sky centered on the 
Praesepe cluster (see Figured]). The question we now pose is whether these are all the RRab 
stars that could have been detected in this distance range given the available PTF data. To 
answer this question we need to know the success rate of our period-finding algorithm as the 
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period estimation is the most critical step in the selection of RR Lyrae stars. If the period 
of pulsation is not correctly estimated, the phased light curve will not be identified as that 
of a RR Lyrae star. 

To estimate the success rate of our period-finding algorithm, we have simulated 100 
timeseries by randomly sampling 133 data points from a template RRab light curve with a 
maximum brightness Rq = 20 mag (Rq of the faintest RRab star in our sample), amplitude 
A R = 0.3 mag, and period of 0.58 days. Random sampling of a template light curve is a 
reasonable choice since the observed light curves are almost uniformly sampled (see Figure H]) . 
The amplitude of Ar = 0.3 mag represents the lower range of amplitudes that RRab stars 



can have (e.g., see Figure 16 in lSesar et al.ll2010af ) and the number of data points corresponds 



to the number of observations of the least observed RRab star in our sample (see Tabled]). 
To make the light curves realistic, we added Gaussian noise to sampled magnitudes using a 
model that describes the mean photometric uncertainty in the PTF i?-band as a function of 
magnitude. We then ran mock time-series through our period-finding algorithm and counted 
time-series where the estimated period was within 1% of the true period. 

Even for this worst-case scenario, we find that the periods are successfully recovered in 
better than 95% of cases and conclude that our selection procedure has recovered almost 
every RRab star in the Praesepe region between 77 and 96 kpc from the Sun; at most one 
RRab star may have been missed. 



5. Spectroscopic Observations 

The group of 8 RR Lyrae stars detected in the Praesepe region at 86 kpc from the Sun, 
and described in the previous section, is remarkable as it is the densest and most numerous 
group that we have detected at these distances in our preliminary search. However unlikely, 
this group could still be a chance association of 8 distant stars and line-of-sight velocities 
are needed to test whether the stars in this spatial group also form a moving group. 



5.1. Instrument Setup and the Observing Strategy 



The spectros copic observations were obtained using two instruments: the Double Spec- 
trograph (DBSP; bke fc Gunnlll982h mounte d on the Paloma r 5.1-m telescope and the Low- 
Resolution Imager and Spectrograph (LRIS; lOke et al.lll995l ) mounted on the Keck I tele- 
scope. A 1" wide slit, a 600 lines mm -1 grating, and a 5600 A dichroic were used with 
both instruments. This instrument setup covers a spectral range from 3800 A to 5700 A 
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and results in resolutions of R = 1360 and R = 1760 for blue channels of DBSP and LRIS, 
respectively. 

The target RR Lyrae stars were observed over the course of several nights in sets of one 
or two consecutive exposures followed by observations of a FeAr (DBSP) or HgCdZn (LRIS) 
arc lamps. In order to avoid the discontinuity in the radial velocity curve near the maximum 
light, the observations were scheduled to target stars between phases of 0.1 and 0.85 of their 
pulsation cycle, with earlier phases being preferred as the stars are brighter then. 

In add i tion, on one night we also targeted a RRab star (AV Peg) previously observed 



by iLaydenl (119941 ) . This star was observed in order to validate the metallicity and systemic 
velocity measuring procedures described below. The log of spectroscopic observations is 
provided in Table [2J 



5.2. Data Reduction and Calibration 

All data were reduced with standard IRAF^I ro utines, and s pectra were extracted us- 



ing an optimal (inverse variance- weighted) method (jHornd Il986l ). The spectra of objects 
observed in two exposures were reduced and calibrated independently and then combined. 
The wavelength calibration of DBSP spectra was done using a single set of 49 FeAr lines, 
while the LRIS spectra were calibrated using a set of 13 HgCdZn lines. The highest root- 
mean-square (rms) scatter of the fit of the wavelength calibration was 0.1 A (~ 6 km s _1 at 
4750 A) for DBSP spectra and 0.05 A (~ 3 km s" 1 at 4750 A) for LRIS spectra. 

To correct for the possible wavelength shift in DBSP spectra during an exposure (e.g., 
due to instrument flexure), we measure the wavelengths (A Q b s ) of three [Hgl] sky lines and 
offset the wavelength zero-point by the average of AA = A Q b s — Ai a b, where Aj a b = 4046.565, 
4358.335, and 5460.750 A are the laboratory wavelengths of [Hgl] sky lines when observed 
in air. The uncertainty in the wavelength zero-point (<J wave in units of A) is estimated using 
the standard deviation of AA. 

To correct LRIS spectra we use the [OI] sky line at Ai a b = 5577.34 A. As this line is 
outside the range of arc lines used to calibrate LRIS spectra (the reddest arc line is at ~ 5460 
A), it is possible its position, and therefore the wavelength zero-point of LRIS spectra, may 
be more uncertain if the shape of the wavelength solution near 5577.34 A chan ges significantly 
throughout the night. This would correspond to a distortion of the LRIS focal plane, not to 
flexure, which is essentially a differential offset along this plane. To estimate the uncertainty 



e http : //iraf .noao . edu 
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in the shape of the wavelength solution near 5577.34 A, we use the following procedure. 
First, we evaluate on a fixed grid all wavelength solutions obtained from arc images taken 
during a night. The solutions are offset to the same wavelength zero-point by subtracting 
the average difference between a solution and some reference solution. Then, we evaluate 
all solutions at a single grid point for which the reference solution returns 5577.34 A. We 
find the rms scatter around this point (5577.34 A) to be very similar to the rms scatter of 
the fit of wavelength solutions (~ 0.05 A or 2-3 km s _1 ). This means that the shape of the 
wavelength solution near 5577.34 A does not vary significantly, even though this region is 
outside the range of arc lines used in wavelength calibration. Therefore, for LRIS spectra 
we adopt the rms scatter of the fit of the wavelength calibration as the uncertainty in the 
wavelength zero-point (a wave ~ 0.05 A). Later in Section |5T3| we will use a wave (in units 
of A) to estimate the zero-point uncertainty in measured radial velocities (a zp t in units of 
km s _1 ). 

The signal-to- noise ratio (S/N) of the final spectra at 4750 A ranged from 10 to 20 for 
DBSP spectra and 20 to 30 for LRIS spectra. Figure [3] shows a few examples of observed 
DBSP spectra, ordered from the highest to lowest S/N. The Balmer lines (starting at H/3) 
and [Ca II] H and K lines are easily seen in the spectra. The CH G-band absorption feature 
at 4300 A is relatively weak in spectra of RR Lyrae stars, indicating that the RR Lyrae stars 
shown in Figure |3] are likely more metal poor than the ELODIE template star HD 22879 
([Fe/H] = -0.84 dex). The ELODIE template star HD 22879 is a F9 dwarf (T eff = 5841 K) 
and is one of the radial velocity standard stars used in the next section. 



5.3. Line-of-Sight Velocities 

Heliocentric radial velocities, v r , were obtained by cross-correlating our spectra with 
spectra of 11 F2 to F 9 dwarf stars (hereaft er template spectra, see Table [3]) taken from 



the ELODIE catalog!! ( IMoultaka et al.ll2004l ). The template spectra were degraded to the 
resolution of our observations and were selected to encompass the range of spectral types 
covered by RRab stars as they pulsate. Expanding the list of template spectra to include 
more stars and more spectral types did not significantly reduce the uncertainty of measured 
radial velocities. 

The cross-correlation was performed using the FXCOR task in IRAF, with template 
and science spectra normalized to the pseudo-continuum prior to cross- correlation. For 
our relatively low resolution spectra, the narrowest cross- correlation peaks were obtained 



'http : //atlas . obs-hp . f r/elodie/ 
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by restricting the correlation region in template and science spectra to cores and wings of 
Balmer lines (spectral regions used in cross-correlation are shown in Figure [3]). Addition of 
other spectral regions or Fourier filtering did not significantly improve the cross-correlation. 

The radial velocities were measured by fitting a Gaussian to the top 25% of the cross- 
correlation peak and the final v r radial velocity of an object was calculated as the average 
of the 11 v r values obtained from 11 cross-correlations with ELODIE templates, where each 
v r value was weighted by its cross-correlation error a cc . The cross-correlation error a cc 
reported by FXCOR includes the uncertainties due to broadening of spectral features during 
an exposure, random errors in wavelength calibration, and errors due to ELODIE template 
mismatch. Therefore, radial velocities from better template fits contribute more towards the 
final radial velocity value. 

In addition to the cross-correlation error, there is an error in radial velocity due to the 
uncertainty in determining the wavelength zero-point. We calculate this zero-point error as 



On 



a 



zpt 



4750 



(7) 



in units of km where a wave is the uncertainty in wavelength zero-point in units of 
Angstrom (see Section 15.2 j) . The final v r radial velocities and their uncertainties (cross- 
correlation and zero-point) are listed in Table [2j 



Following lHawley fc Barnes! ( 119851 ). the systemic (center-of-mass) velocity of RR Lyrae 
stars, Vj, was obtained by fitting the radial velocity curve of the RRab star X Ari to observed 
heliocentric velocities. We used Layde n's (1994) parametrization (see his Figure 8) of the 
radial velocity curve that lOkd ( 119661 ) measured from the H7 line. This template is an 
appropriate choice for fitting radial velocities measured using Balmer lines (i.e., as done 
in this work). Radial velocity t emplates based on measurements of metallic lines, such as 
the one presented in iLiul (1199ll ). would not be appropriate here a s radial ve locity curves 
of metallic and Balmer lines have different shapes and amplitudes (jOkdll966l ). The X Ari 
template was shifted in velocity to match the measurement for each star at the corresponding 
phase, and the systemic velocity is that of the shifted template at phase 0.5. Two examples 
of template fits are shown in Figure HI 

Finally, the error calculation for the systemic velocity i ncludes terms tha t take into 
account errors in the model radial velocity curve. Following I Vivas et al.l ( 120051 ) (see their 
Section 3), we use the expresion given below in which the third term accounts for likely 
uncertainties in the phase where the velocity curve passes through the systemic velocity, and 
the fourth term is related to possible variations in the slope of the velocity curve: 



= + °L + (H9.5 x 0.1) 2 + (23.9A# 
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where A<p = <p b s — 0.5, and 4> b s is the phase of observation. The a cc and a zpt are the cross- 
correlation and zero-point error, respectively. The contributions of model uncertainties are 
of the order of ~ 15 km s _1 and are higher than the observed radial velocity errors [a cc and 
cr zpt added in quadrature). Therefore, to reduce the uncertainty in the systemic velocity we 
need to reduce the uncertainty introduced by the model, that is, we should observe the radial 
velocity curve more than once. 

When there is more than one v 7 estimate available for a RRab star due to multiple 
measurements, we simply calculate the weighted average of t> 7 , (t> 7 ), and its associated 
uncertainty. Finally, we correct the heliocentric systemic velocities (v 7 ) for the solar motion 
with respect to the Galactic center using 

(vgsr) = (%} + v u cos b cos I + {vy + vlsr) cos b sin / + v w sin 6, (9) 



wher e (vu, ty, v w) = (10.0,5.2,7.2) km s 1 and vlsr = 220 km s 1 (IBinney k, Merrifield 



19981 ). where vlsr is the velocity of the local standard of rest. These velocities in the 
Galactic rest frame {(v gsr )) are adopted in the rest of this paper. 

To validate the above procedure, we measured the radial velocity of AV Peg . The light 



curve parameters for this star were taken from the GEOS RR Lyrae database (ILe Borgne 



20071 ). After correcting the observed heliocentric velocity for pulsation using the model 



radial velocity curve, we obtained (t> 7 ) = —64.8 ± 18.7 km s \ Thi s value agrees within the 



uncertainties with the (v 7 ) = —58 ± 1 km s 1 estimated by iLaydenl ( 19941 ) (see his Table 9) 



5.4. Metallicities 



Spectroscopic metallicities were measured following the method and calibration of |L ay den 



( 1l994j ) which involves comparing the pseudo-equivalent width of [Call] K line, W(K), against 



the mean pseudo-equivalent widths of (3, 7, and 5 Balmer lines, W(H). 

The pseudo-equivalent widths of the [Ca II] K line and Balmer lines were measured using 
the EWIMH program^ (written by A. Layden) from DBSP and LRIS spectra normalized to 
the pseudo-continuum. Measured pseudo-equivalent widths (W'(K), W'(H5), W'{H^), and 
W'(H(3)) were transformed to the Layden's (1994) pseudo-equivalent width system using 

W(K) = 1.23W'(K) - 0.44 (10) 
W(H5) = 0.84W'(H5) + 1.64 (11) 



s http : //phy sics . bgsu . edu/ "layd en/ ASTR0/DATA/EXP0RT/EWIMH/ ewimh . htm 
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W(H^) = 1.13W'(Hj) - 0.63 
W(HP) = l.Q2W'(Hp) + 0.98 

for pseudo-equivalent widths measured from DBSP spectra, and using 

W(K) = 1.16W'(K) - 0.39 
W{H8) = l.02W'(H8) + 1.39 
W(Hi) = 1.07W'(Hj) - 0.45 
W(H0) = 1.01W'(Hf3) + 1.18 



(12) 
(13) 

(14) 
(15) 
(16) 
(17) 



for pseudo-equivalent widths measured from LRIS sp ectra. To der ive these equations we 
observed 5 of the equivalent width standard stars of iLaydenl (119941 ) with DBSP and four 
more with LRIS; se e Table H I Com paring our measurements of their pseudo-equivalent widths 
with those listed in ILaydenl (119941 ) yielded Equations [TO] to [TT1 Af ter the transformation, the 
W(K) was corrected for interstellar [Call] absorption using the iBeerd (119901 ) model. This 
correction decreases the measured metallicities by ~ 0.15 dex. 

The spectroscopic metallicity was calculated as 

= W(K)-a-W { H) 
L ' 1 c + dW{H) ' V ; 

where a = 13.858, b = -1.185, c = 4.228, d = -0.3 2 (see Table 8 in iLavdenl fl99^ . 
Equation IT81 was obtained by inverting Equation 7 from ILaydenl (I1994I ). The metallicities 
obtained using this method are listed in Table [2j 

For the most part, [Fe/H] estimates obtained from multiple observations of the same 
star are consistent within 0.1 to 0.2 dex, even for RR1 which was observed with different 
instruments (DBSP and LRIS). This scatter is con sistent with the systematic uncertainty 
of this method, which is ~ 0.15 dex ( ILaydenl Il994l ). The exception to this is star RR3 for 
which we measured [Fe/H] = — 1.5 dex from a DBSP spectrum and [Fe/H] = —2.0 dex from 
a LRIS spectrum. If we assume [Fe/H] = — 1.8 dex as the metallicity of this star (the average 
of -1.5 and -2.0), we find that the two measurements are still within 2a of the average value, 
which is statistically not improbable. While it may be tempting to assign lower significance 
to DBSP metallicity measurements due to lower S/N of DBSP spectra (10 to 20 vs. > 20 
for LRIS spectra), multiple metallicity measurements of RR5 argue against that. RR5 has 
been observed with DBSP over three nights, its spectra have S/N in the 10 to 20 range, and 
yet the standard deviation of its metallicity measurements is 0.1 dex, consistent with the 
systematic uncertainty. 

To validate the zeropoint o f our metallicity scale, we measured [Fe /H] of RRab st ar AV 
Peg observed by ILaydenl ( J1994J ) . We obtained [Fe/H] = 0.1 dex while ILaydenl (119941 ) cites 
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[Fe/H] = —0.1 dex. These two values are roughly within the systematic uncertaint y of this 
metho d, indicating t hat our metallicitie s are on the same zeropoint syste m as the lLayden 
(119941 ) stars (i.e., the IZinn fc Westlll984j globular cluster abundance scale; lLaydenlll994[ ). 



6. Discussion 

In this section we address the following questions. 

1. What is the mean velocity and velocity dispersion of our sample, and how do these 
values compare to the literature? 

2. Is the observed distribution of velocities a Gaussian or a non-Gaussian distribution? 

3. Are the observed velocity groups bound or unbound systems? 

4. What does the metallicity of their stars say about the progenitor (s)? 

5. Are the observed velocity groups related to known halo substructures? 

The galactocentric rest-frame velocities (v gsr ) listed in Table |2]indicate that the observed 
RR Lyrae stars can be split into two velocity groups, one moving at v gsr = 78.0 ±5.6 km s _1 
and the other one at v gsr = 16.3±7.1 km s~~ , where v gsr is the weighted mean of (v gsr ) values 
listed in Table HJ The uncertainty in v gsr is calculated as the standard error of the weighted 
mean. We tentatively name these moving groups Cancer group A and B. The weighted mean 
velocity for the entire sample is v gsr = 53.5 ± 4.4 km s~ . 

To estimate the velocity dispersion of these moving groups, we use the following Monte 
Carlo procedure. We generate 1000 sets where each set contains four mock line-of-sight ve- 
locities, v moc k (four values because there are four RRab stars in each moving group). A v moc k 
value is generated by drawing a single value from a Gaussian distribution centered on the 
(v gs7 ) value of a star (e.g., 75.6 km s _1 for RR1) and having a width equal to the uncertainty 
on (v gsr ) (e.g., 8.5 km s" 1 for RR1). The weighted standard deviation is calculated and 
stored for each set, with v moc k values weighted by uncertainty on corresponding (v gsr ) values 
listed in Table [21 The velocity dispersion is then the average value of the distribution of 
1000 weighted standard deviations and the uncertainty of velocity dispersion is the standard 
deviation of this distribution. Using this procedure we find that the two moving groups have 
velocity dispersions of cr Vgsr = VIA ± 5.0 km s _1 (group A) and <J Vgsr = 14.9 ± 6.2 km s -1 
(group B). The velocity dispersion of the entire sample is <J Vgsr = 37.7 ± 5.0 km s -1 . 
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The measured velocity dis persion for th e full s ample is inconsistent with the halo velocity 
dispersion profile obtained by iBrown et al.l (120101) usin g a sample of 910 distant halo stars. 



At 92 kpc, the extrapolation of the IBrown et all (120101 ) model (see their Equation 6) would 
describe the velocity distribution of stars as a 68 km s _1 -wide Gaussian centered on km s _1 . 
The mean velocity of our full sample is at least 10<r away from km s _1 , and the velocity 
dispersio n of our full sample is at least 5<r lower than the velocity dispersion extrapolated 
from the IBrown et al.l (120101 ) model. One possible explanation of this discrepancy is that 
Brown et al.l (120101 ) Equation 6 cannot be extrapolated outside the 16 < Rqc < 64 kpc 
range within which the relation is assumed to be valid (IBrown et al.ll2010l ). A more likely 
explanation is that we are observing two moving groups with their own internal kinematics 
(mean velocity and velocity dispersion). If that is the case, a comparison of the predicted 
and measured velocity dispersion for the full sample is not a sensible thing to do. 

The question of whether the distribution of velocities listed in Table [2] is a Gaussian or 
a non-Gaussian distribution is an important one, as deviations from normal ity are usually 



interpreted as a signature of velocity group s in the smooth halo component (IHarding et al. 



2001; Duffau et al. 



2006 



Vivas et al. 



20081 ). T o test for the pre s ence of velocity groups in 



the halo, IHarding et al.l (120011 ) recommend the IShapiro fc Wilkl (119651 . SW) statistical test 
of normality to be applied to velocity histograms. This test is sensitive to many different 
deviations from the Gaussian s hape and does not de pend on the choi ce of mean or dispersion 



of th e normal distribution (see IHarding et al.ll200ll Section 5.1 and iD'Agostino &: Stephens 



1986 for more details on the test) 



Following the procedure detailed in ISesar et al.l (j2010bl . see their Section 3) that takes 
into account uncertainties in velocities, we use the SW test and find that there is a 37% 
probability that the observed sample of velocities was drawn from a single Gaussian distri- 
bution. However, the SW test does not have the ability to assess the rather peculiar nature 
of the shape of the observed distribution (i.e., it cannot distinguish the particular nature of 
non-Gaussianity). To include this information, we consider the 37% chance that the mea- 
sured velocities were drawn from a Gaussian distribution, and estimate the probability that 
this chance would produce the observed distribution seen in Figure [51 We generate 100,000 
samples and for each sample draw 8 velocities from a Gaussian centered on 53.5 km s _1 
and 37.7 km s _1 wide (i.e., the velocity distribution of the full sample of RR Lyrae stars 
observed in the Cancer region). Each sample is then divided into stars with velocities greater 
and smaller than 53.5 km s _1 ("positive" and "negative" groups), and the mean velocities 



and velocity dispersions are calculated for each group 



'J pos i ^neg > 



@pos J 

and a neg ). In only 



~ 1.6% of generated samples, we find that the "positive" and "negative" groups have mean 
velocities and velocity dispersions in ranges covered by Cancer groups (72 < v pos < 84 and 



7 < a pos < It 



and 9 < v neg < 24 and 8 < a neg < 22). The percentage of generated 
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samples with mean velocities and velocity dispersions in ranges covered by Cancer groups 
is even lower (0.02%), if the samples are drawn from a 68 km s _1 -wide Gaussian centered 
on km s _1 ( i.e., the velocity distribution at 92 kpc as extrapolated from Equation 6 by 



Brown et al.ll2010l ). Therefore, even if the observed velocities were drawn from a Gaussian 
distribution (a 37% chance), it is unlikely they would create the distribution seen in Figure El 
as the chance of that happening is only about 0.6% (0.37 x 0.016 ~ 0.006), if the velocity 
distribution at 92 kpc is a 37.7 km s _1 -wide Gaussian centered on 53.5 km s" 1 , or about 
0.07% (0.37 x 0.000 2 ~ 7 • 10~ 5 ), if the velocity distribution follows the extrapolation of the 
Brown et al.l (120101 ) model. The next simplest explanation for the observed distribution of 



velocities is that we are observing two velocity groups, moving at v 



si-r- 



and v 
and <7 V 



gsr 



16.3 ±7.1 km s 1 and having velocity dispersions of a v 
15.0 ±6.1 km s~ x . 



78.1 ± 5.7 km s" 
12.4 ± 5.0 km s" 



These estimate d velocity dispersio ns are consistent with velocity di spersions of dwarf 
spheroidal galaxies (IWalker et al.ll2007l . e.g., Leo II) and globular clusters (lHarrislll996l . 2010 
edition, e.g., NGC 362). Due to significant uncertainities in the velocities, these velocity 
dispersions are most likely upper limits on the intrinsic velocity dispersions of the Cancer 
groups. While the velocity dispersions of the two groups are consistent with velocity dis- 
persions of dwarf spheroidal galaxies and globular clusters, the groups are spatially quite 
extended (the maximum separation between members is several kpc) making it very un- 
likely that they are bound systems, and are more likely to be debris of tidally disrupted 
dwarf galaxies or globular clusters. 

The metallicity of these groups may provide more information about their progenitors. 
The first group has a me dian metallici t y of [ Fe/H] = — 1.6 dex and the second group has 
[Fe/H] = —2.1 dex on the IZinn fc Westl (119841 ) globular cluster abundance scale, indicating 



that the progenitors of these two substructures are systems with old, metal-poor populations. 
The metallicity spread (standard deviation) in both groups is rather large (~ 0.4 dex) and a 
few times greater than the uncertainty in individual metallicity estimates (~ 0.15 dex). Such 
a large dispersion in [Fe/H] indicates self-enrichment and implies the progenitor structures 
were dwarf galaxies and most likely not globular clusters. Using the metallicity-luminosity 
relation by iKirby et al.l ( 1201 ll ) (see their Equation 9), we can place a lower limit on the 
luminosity of group A and B progenitors. Taking into account the intrinsic spread in this 
relation, we obtain L > 5.5 x 10 5 L for the group A progenitor and L > 9.1 x 10 3 L Q for 
the group B progenitor. The luminosity limit of L > 5.5 x 10 5 L indicates that the group 
A progenitor was likely more lumin ous than one of th e classical dwarf spheroidal galaxies, 
such as Leo II (e.g., see Figure 3 by lKirby et al.l 1201 ll ). 
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However, we do note that there are only four stars in each group and that the metallicity 
spread can change drastically if a single measurement is not considered when estimating the 
spread. For example, if the metallicity of RR8 is not taken into account when estimating 
the metallicity spread of the second group, the new spread would be ~ 0.1 dex which would 
be consistent with measurement uncertainties and consistent with a single stellar population 
(i.e., a globular cluster as the progenitor). Additional observations of RR Lyrae stars on the 
metal-poor and metal-rich ends (e.g., RR2, RR3, and RR8) would be very useful as they 
would further test the validity of the estimated metallicity dispersions. 



Sharma et al.l (120101 ) have found two overdensities of M giants extending into this region 
of the sky, one at 22.6 ± 11 kpc (their group A13) and the other one at 96.6 ± 48 kpc (their 
group A9). Both of these groups have density peaks (centers) that are offset by at least 20° 
from the location of Cancer groups discussed here, with the group A13 peaking at 1 ~ 144° 
and b ~ 31° and gr oup A9 peaking at 1 y 188° and b ~ 20°. The Cancer groups are likely not 
associated w ith the ISharma et al.l (120101 ) group A13 based on estimated distances of M giants 
provided by ISharma et al.l (120101 ) (22.6 ± 11 kpc, see their Table 2). Based on the distance 
alo ne (96.6 ± 48 kpc), i t is not impossible that the Cancer groups may be associated with 
the ISharma et al.l (120101 ) group A9. Ho wever, we do n o te tha t the Cancer groups discussed 
in this work are on the very edge of the ISharma et al.l (120101 ) group A9 (see their Figure 7) 
and that the distance of group A9 has uncertainty of ~ 50 kpc. 



Newberg et al.l (120031 ) have found an overdensity of A-type stars about 10° north in 
declination from Cancer groups and at a similar distan ce (~ 90 kpc from the Sun; see their 
Figure 1). At a similar location. iRuhland et al.l (120111 ) have found "a very faint indication 
of an overdensity" of blue horizontal branch (BHB) stars. Both studies speculated that 
the overdensities they observed were associate d with the Sagitta r ius ti dal debris. However, 
recent modelling of the Sagittarius streams by lLaw fc Majewskil (120101 ) suggests that there 
should not be any debris at these distances. It is unclear whether these overdensities and 
C ancer groups ar e relat ed d ue to lack of sp e ctrosc opic follow-up of BHB and A-type stars 



in 



Newberg et al.l (120031 ) and lRuhland et al.l (120111 ) overdensities 



The surface density of RRab stars in Cancer groups A and B seems to be similar to that 
of the Pisces Overdensity /Stream (~ 0.3 RRab per d eg 2 ), a halo subs t ructure located at a 
similar distance (~ 85 kpc) from the Galactic Center (jSesar et al.ll2007l ; IWatkins et al.ll2009 



Sesar et al.ll2010al ). In addition, both of these substructures seem to have tw o velocity groups 



(the Pisces Overdensity has groups moving at 50 km s 1 and —52 km s ; iKollmeier et al. 



20091 ; ISesar et al.l l2010bl ) . The morpholo gy of both substructures i s reminiscent of debri s 



clouds observed in some simulat ions (e.g.. iBullock fc Johnston! 120051 ; iJohnston et al.l 120081 ) . 



According to a recent study by Johnston et al. ( 20121 ). the clouds represent debris that 



is 



-17- 



slowly turning around at apocenters of orbits that are typically more eccentric than orbits of 
debris streams, such as the Sagittarius stream. If so, then "there must be low-density stellar 
streams movi ng; between these apo centric clouds and passing through the inner Galaxy at 
high speed." ( IJohnston et al.ll2012[ ). By following up RR Lyrae stars selected from PTF and 
other surveys, it will be possible to identify such low density streams, if they exist, and to 
explore their relationship with Cancer groups A and B. 



7. Conclusions 



We confirm the existence of two kinematic groups in the direction of the Cancer con- 
stellation (R.A ~ 129° and Dec ~ 20°, or 1 ~ 205° and b ~ 32°), located at 92 kpc from 
the Galactic center (86 kpc from the Sun). These groups, tentatively named Cancer groups 
A and B, are moving at v^ sr = 78.0 ± 5.7 km s" 1 (Cancer group A) and Vg Sr = 16.3 ± 3.8 
km s _1 (Cancer group B). The groups have velocity dispersions smaller than 15 km s _1 , are 
spatially extended (about several kpc), metal-poor (median metallicities of [Fe/H] A = —1.6 
dex and [Fe/H] B = —2.1 dex), and have a metallicity spread of ~ 0.4 dex. These results 
suggest that the observed groups are debris of tidally disrupted dwarf galaxies, possibly near 
the apocenters of their orbits. Whether these groups are related to known substructures in 
the Galactic halo is unclear at this point, and answering this quest ion may require exten sive 
orbit modeling an d comparisons with simulations (e.g., as done by lJohnston et al.ll2012l and 
Carlin et allboij l. 



Observations of two RR Lyrae stars (RR9 and RR10 in Tables [T] and [2]) obtained after 
this paper was submitted for peer-review, may help with the modeling of Cancer groups' 
orbits. The two RR Lyrae stars have velocities of 88 km s -1 and 38 km s _1 that are consistent 
(to within the uncertainties) with the mean velocities of Cancer groups A and B, and are 
at similar distances (~ 76 kpc from the Sun). They are offset ~ 25° east of the Cancer 
groups and may indicate an eastward extension of these groups. However, due to the sparse 
coverage in PTF of the sky between Cancer groups and RR Lyrae stars RR9 and RR10, we 
are unable to verify at this point whether RR9 and RR10 are truly related to the Cancer 
groups or not. 

Initially, this work was motivated by a hypothesis that distant halo substructures may 
be found by simply looking for distant spat i al gro ups o f RR Lyrae stars. Based on this 
work and previous work by iKollmeier et all (j2009[ ) and ISesar et all (l2010b[ ). we conclude 
that this is indeed an efficient approach to finding and following-up halo substructures. So 
far, all of the distant (galactocentric distances greater than 80 kpc) RR Lyrae stars located 
in spatial groups have proven to be members of a halo substructure (e.g., Cancer groups in 
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this w ork and RR Lyrae stars in the Pisces Overdensity; iKollmeier et al.ll2009l ; ISesar et al. 



2010 bf). Since the density profile of th e relatively smooth, inner halo steepens beyond 30 kpc 



( ISesar et all 1201 It iDeason et al.ll201ll ). it may be that the majority, if not all, of RR Lyrae 
stars beyond ~ 30-40 kpc are part of some halo substructure. Spectroscopic followup of 
distant RR Lyrae stars not studied in this work may provide more data to support or refute 
this hypothesis, and we plan to follow this strategy with other RR Lyrae stars observed by 
the Palomar Transient Factory. 
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Table 1. Positions and Light Curve Parameters of RR Lyrae Targets 



ID" 


IAU name 


R.A. b 


Dec 6 


distance c 


A R 


R o 


Period^ 


HJD» 


N obs 






(deg) 


(deg) 


(kpc) 


(mag) 


(mag) 


(d) 


(d) 


RR1 


PTF1J084844.61+202430.8 


132.185880 


20.408611 


87.1/93.2 


0.887 


19.84 


0.515257 


55303.686299 


291 


RR2 


PTF1J084313. 05+193901.1 


130.804370 


19.650306 


78.2/84.4 


0.889 


19.81 


0.556328 


55271.737376 


414 


RR3 


PTF1J083851. 75+221641. 5 


129.715656 


22.278200 


85.9/92.2 


0.766 


19.82 


0.555729 


55590.910809 


160 


RR4 


PTF1J082914.96+185607.4 


127.312372 


18.935392 


86.0/92.4 


0.629 


20.01 


0.594511 


55667.607592 


193 


RR5 


PTF1J083636.44+200954.9 


129.151846 


20.165259 


96.0/102.3 


0.855 


19.92 


0.610057 


55594.666378 


406 


RR6 


PTF1J083339.04+192641.5 


128.412679 


19.444856 


92.0/98.4 


0.697 


19.97 


0.633480 


55823.965992 


271 


RR7 


PTF1J083033. 30+220217.9 


127.638788 


22.038294 


77.0/83.5 


0.703 


19.64 


0.545361 


55594.799336 


400 


RR8 


PTF1J082333. 11+213414.6 


125.887964 


21.570710 


84.2/90.8 


0.912 


19.96 


0.517656 


55542.999049 


133 


RR9 


PTF1J103051. 69+202330.8 


157.715375 


20.391878 


75.5/79.3 


0.914 


19.35 


0.544278 


55515.539757 


82 


RR10 


PTF1J102040. 66+212315.1 


155.169417 


21.387514 


76.1/80.2 


0.759 


19.52 


0.611088 


55244.324620 


96 


AV Peg 


AV Peg 


328.011644 


22.574827 


0.7/7.9 


1.040 


n.a. 


0.390381 


55833.513400 


n.a. 



Note. — The first 8 RR Lyrae stars are located in the Cancer region, the next two are offset by ~ 25° east from the Cancer region and 
suggest an eastward extension of Cancer groups, and AV Peg is a known RR Lyrae stars used for validation of various procedures used in 
this work. 

a Short name. 

b Equatorial J2000.0 right ascension and declination from SDSS DR8 catalog. 
c Helioccntric/galactocentric distance calculated using spectroscopic mctallicity. 
d Light curve amplitude in the Mould- R band (V-band amplitude for AV Peg). 
°i?-band magnitude at the epoch of maximum brightness. 
Period of pulsation. 

s Rcduccd Heliocentric Julian Date of maximum brightness (HJD - 2400000). 
h Number of observations in the Mould- R band. 



Table 2. Observing Log, Velocities, and Metallicities 



ID" HJD s 6 pectrum Exposures- Phase 1 * v T ± a cc ± a e zpt v- i ±a S 1 (v^ (v gsr ) h [Fe/H]' 





I'd"! 


(s) 




(km s ) 


(km s ) 


(km s ~) 


(km s ) 


(dcx) 


RR1 


55866.028827 


2 X 1200 


0.38 


151.2 ±8.6 ±4.0 


165.5 ± 15.5 


159.9 ±8.5 


75.6 ±8.5 


-1.5 


RR1 


55926.934245 


1 x 1800* 


0.59 


170.0 ±5.7 ±2.5 


159.2 ± 13.6 






-1.7 


RR1 


55952.995976 


1 x 1800* 


0.17 


115.8 ±4.8 ±3.2 


155.2 ± 15.4 






-1.7 


RR2 


55865.996308 


2 x 1200 


0.18 


117.8 ±6.8 ±5.7 


156.0 ± 16.7 


156.0 ± 16.7 


69.1 ± 16.7 


-0.9 


RR3 


55917.913149 


2 x 1800 


0.42 


155.0 ±4.2 ±6.9 


164.6 ± 14.6 


157.0 ±9.8 


79.9 ±9.8 


-1.5 


RR3 


55953.024432 


1 x 1800* 


0.60 


162.8 ±4.0 ±3.2 


150.9 ±13.2 






-2.0 


RR4 


55916.996225 


2 x 1800 


0.49 


181.5 ±9.8 ±10.9 


182.7 ±18.9 


182.7 ± 18.9 


93.9 ± 18.9 


-1.4 


RR5 


55895.934069 


2 x 1200 


0.84 


129.5 ± 7.1 ± 16.8 


98.2 ±23.3 


91.3 ± 12.4 


6.5 ± 12.4 


-2.6 


RR5 


55922.944423 


2 x 1800 


0.11 


47.9 ±4.9 ±7.7 


94.5 ± 17.7 






-2.4 


RR5 


55926.018096 


2 x 1800 


0.15 


34.2 ±7.5 ±19.9 


76.0 ±25.8 






-2.5 


RR6 


55926.967532 


1 x 1800* 


0.60 


117.3 ±5.4 ±2.5 


105.4 ±13.6 


105.4 ± 13.6 


18.1 ± 13.6 


-2.1 


RR7 


55917.864566 


2 x 1800 


0.39 


87.5 ±3.9 ±9.9 


100.7 ±16.2 


100.7 ± 16.2 


23.0 ± 16.2 


-2.2 


RR8 


55888.997603 


2 x 1800 


0.39 


88.1 ±8.2 ±2.2 


101.3 ±14.9 


101.3 ± 14.9 


22.4 ± 14.9 


-1.2 


RR9 


56037.732088 


2 x 1800 


0.50 


109.8 ±16.8 ±3.2 


109.8 ± 17.1 


109.8 ± 20.9 


38.0 ± 20.9 


-2.0 


RR10 


56037.175892 


2 x 1800 


0.44 


152.1 ± 11.8 ± 5.7 


159.3 ± 13.1 


159.3 ± 17.8 


88.6 ± 17.8 


-1.9 


AV Peg 


55917.604262 


1 x 270 


0.41 


-75.5 ±7.8 ±11.9 


-64.8 ±14.2 


-64.8 ± 14.2 


135.1 ± 14.2 


0.1 



a Short name. 

b Reduced Heliocentric Julian Date when the spectrum was taken (HJD - 2400000). 
-Exposure times (* denotes Kcck-I/LRIS setup). 
d Pulsation phase when the spectrum was taken. 

-Radial velocity (not corrected for pulsations) and its uncertainties (cross-correlation and zero-point). 

f Systemic (center-of-mass) velocity (corrected for pulsations) and its uncertainties (cross-correlation, zero-point, and model errors 
added in quadrature). 

g Weighted average of multiple heliocentric systemic velocities and its uncertainty. 

h Galactocentric rest-frame velocity and its uncertainty. 

'Metallicity. 



Table 3. ELODIE Template Stars 



Name 


ELODIE ID Number 


Spectral Type 


HD 49933 


00282 


F2V 


HD 140283 


00480 


sdF3 


HD 102870 


00414 


F8V 


HD 693 


00004 


F5V 


HD 7476 


00690 


F5V 


HD 13555 


00998 


F5V 


HD 3268 


00010 


F7V 


HD 222368 


00855 


F7V 


HD 19994 


00084 


F8V 


HD 22484 


00092 


F9V 


HD 22879 


00096 


F9V 



Table 4. Equivalenth-width Standard Stars 



Name 


W(K) 

vv V^V 


W'(H<5) 


W'(K-y) 

VV J J 


W'fH/3) 




mean ± sd 


mean ± sd 


mean ± sd 


mean ± sd 


Kopff 27 


1.72 ± 0.06 


10.5 ±0.4 


10.2 ± 0.3 


8.84 ±0.05 


HD 180482 


2.42 ± 0.07 


10.26 ± 0.09 


9.2 ±0.1 


8.7±0.1 


BD+25 1981 


2.79 ± 0.05 


5.0 ±0.2 


5.4 ±0.3 


5.12 ±0.03 


HD 155967 


6.6 ±0.2 


3.0 ±0.3 


3.89 ±0.08 


3.95 ±0.09 


HD 112299 


8.4 ±0.2 


1.9 ±0.2 


2.10 ±0.07 


3.02 ±0.06 


Kopff 27 


1.89 


9.68 


10.44 


8.89 


BD+25 1981 


2.90 


5.47 


6.05 


5.08 


HD 155967 


6.90 


2.75 


3.84 


4.05 


HD 112299 


8.97 


1.64 


2.10 


2.98 



Note. — The horizontal line separates equivalent-width standard stars 
observed with DBSP (over three nights) and LRIS (over a single night). 
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Fig. 1. — Spatial distribution of 8 a&-type RR Lyrae stars studied in this work, which are 
located at distances of 77-96 kpc from the Sun (84-102 kpc from the Galactic Center). Their 
positions and light-curve parameters are listed in Table CD The symbols indicate the halo 
velocity group membership for each of these stars in either the v gsr ~ 80 km s _1 group (blue 
circles) or the v gsr ~ 15 km s~ l group (red squares, see Section [6]). The lines simply connect 
labels and symbols. According to the completeness analysis presented in Section HI these are 
the only RRab stars in this region of the Galactic halo. 
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Fig. 2. — The phased i?-band light curves of three RRab stars discussed in this work (from 
top to bottom: RR5, RR8, and RR1). The light curves are offset for clarity (offsets of 0, 
1.5, and 3 mag from top to bottom ) and the solid lines show best-fit SDSS r-band templates 
constructed by ISesar et al. 
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Fig. 3. — The DBSP spectra of some RRab stars discussed in this work, ordered from lowest 
to highest S/N (top to bottom). An ELODIE template spectrum is shown at the bottom 
of the panel. The DBSP spectra are in the observer's frame and the ELODIE spectrum 
is in the rest frame. All spectra have been normalized to the pseudo-continuum and have 
been offset in flux for clarity. The dashed lines show the rest frame positions of, from left 
to right, [Call] K & H lines and Balmer lines (H<5, H/y, H/3 ). The boxes bracketing Balmer 
lines show the spectral regions used for cross-correlation. The absorption feature redward 
of the H7 line is the [Hg I] 4358 A sky-line that has been oversubtracted during extraction 
of spectra. The dotted line shows the position of the CH G-band at 4300 A. Note how this 
absorption feature is more pronounced in the spectrum of the ELODIE star than in spectra 
of RR Lyrae stars. This indicates that the observed RR Lyrae stars are likely more metal 
poor than the ELODIE star (a F9 dwarf star, with [Fe/H] = -0.84 dex and T eff = 5841 K). 
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Fig. 4.— Fits of the radial velocity template of X Ari (solid line) to two stars in our 
sample. The open circles show individual v r measurements and their error bars indicate the 
uncertainty (cross-correlation and zero-point errors added in quadrature). The red diamond 
symbol at phase 0.5 shows the systemic velocity calculated as the weighted average of v 1 
values, where the v 1 values were obtained by fitting the model template to each of the 
observed radial velocities. The uncertainty in the systemic velocity includes the observational 
errors (cross-correlation and zero-point), and model errors. 
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Fig. 5.— Histogram of (v gsr ) of RR Lyrae stars in Cancer groups. The bin size is 15 
km s -1 , which is about the average velocity error in (v gsr ). The distribution of observed 
velocities seems to be bimodal, with two velocity peaks centered on (v gsr ) = 78 km s -1 and 
(vgsr) = 16 km s _1 , respectively. See Section [6] for a discussion of statistical significance of 
this bimodality. 



